CLAIMS 



We claim: 

1. An iterative method of normalized waveform inversion to obtain a model parameter 
describing one or more physical properties of a medium, the method comprising the 
steps of: 

a) inputting a time domain measurement data set Z)^ (f), = NG; / = 1 NS) , 

where d denotes a data tensor, NG is the number of receivers and NS is the 
number of sources; 

b) means for minimizing a model parameter m below an error bound in the process 
of creating a normalized waveform inversion, wherein said minimizing step 
further comprises: 

i) Fourier transforming the time domain measurement data set Dt(t) to create a 
measurement spectral data set D^, (co) ; 

ii) normalizing the measurement spectral data set D^.(ty) to create a normalized 
data wavefield T^^ico); 

iii) modeling a medium by iterating the model parameter m describing one or 
more physical properties of the medium; by 

(1) minimizing a weighted error, between the normalized data wavefield 
Tj.(q)) and the normalized modeled wavefield T^7(^) th^ response of 

the medium, to a level below the error bound; 

(2) outputting the iterated model parameter m corresponding to the weighted 
error below the error bound, as a minimized model parameter. 

2. In the modeling step of claim 1, the physical property is selected from the group 
comprising: a velocity of sound propagation, density, permeability, porosity, 
resistivity, permittivity, Young's modulus, a component of the tensorial stress-strain 
properties, shear modulus, amount of water present, and amount of oil present. 
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3. In the modeling step of claim 1, the error bound divided by an initial error Eq is 

selected from the group consisting of: below 10"^, below 10'^ , below 10'^ , and below 
10^ 

4. An iterative method of using normalized waveform inversion to obtain a model 
parameter describing one or more physical properties of a medium, the method 
comprising the steps of: 

a) inputting one or more each of NS source and NG measurement spatial locations; 

b) measuring time domain data at each of the NG measurement locations resulting 
from an input waveform at one of the NS source locations propagating through a 
medium, 

i) for each of the NS source locations, 

ii) thereby forming a time domain measurement data set 

D^ji (7 =1'" i = 1 - NS) , where d denotes a data tensor; 

c) Fourier transforming the time domain measurement data set Z>^ (0 to create a 
measurement spectral data set D^,.(6y) having frequency and amplitude 
information for each of the NG measurement locations; 

d) normalizing the measurement spectral data set D'^jfio)) to create a normalized data 

wavefield T^^ico); 

e) modeling the medium using an iterated model parameter m describing one or 
more physical properties of the medium, 

(1) the NS source and the NG measurement spatial locations used as 
respective model input and model response spatial locations contained 
within the model of the medium, 

(2) and initializing the iterated model parameter m with corresponding one or 
more known bulk properties of the medium being modeled, 

ii) said modeHng step comprising: 

iii) creating a measurement model by: 

(1) applying a delta function source collocated with the NS source. 
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(2) modeling the response at the NG measurement locations, using a time 
domain modeling method, to create a synthetic medium response at the 
receiver due to the source, {t),m denoting the model response, 

(3) repeating the applying and modeling steps at each of the NS source 
locations and NG measurement locations until the measurement model is 
full, and 

iv) Fourier transforming the model response (0 » obtain a frequency 
domain synthetic response PJ(^y) ; 

v) forming a normalized modeled wavefield using the frequency domain 
synthetic response T;(ty) = P;(ty)[P,7((y)]"' ; and 

f) minimizing a weighted error, between the normalized data wavefield TjJ(£y) and 
the normalized modeled wavefield THico) of the response of the medium, tp a 
level below an error bound, 

i) said weighted error met by using the iterated model parameter m , known as 
the a minimized model parameter m , 

5. The method of claim 4 wherein the modeled normalized wavefield T^^ict)) is 

independent of the input waveform applied to any of the NS source locations. 

6. The method of claim 4 further comprising the step of: 

a) storing the minimized model parameter w in a computer-readable medium. 

7. The method of claim 4 further comprising the steps of: 

a) displaying the minimized model parameter w as a graphical representation of 
one or more of the physical properties of the medium measured, 

b) storing an output image represented by the minimized model parameter m in a 
computer-readable medium. 

8. The method of claim 4 further comprising the step of; 

a) storing as a computer program in at least one computer-readable medium, 

b) said steps of claim 4. 
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9. The method of claim 4 wherein said weighted error results from a weighing matrix 

that is an inverse of the standard deviation of the normalized data wavefield 

10. The method of normalized waveform inversion of claim 9 wherein said minimizing 
step weighted error (p{m) is calculated by a root mean square of the weighing matrix 

multiplied by the difference between the normalized data wavefield T^i((o) and 

the normalized modeled wavefield X; (o)) , or <Z> (m ) = || (x;* (o)) - X^^ (^y))||^ 

sunmied over all frequencies, / sources, and j receivers. 

11. The method of claim 4 wherein said normalization step is further comprised of the 
steps of: 

a) choosing one of the measurement spectral data set D^,.(ty) at one of the 
measurement locations, such as D^. [o)) , 

(1) said chosen measurement spectral data set DJ^. (co) known as a reference 
measurement; 

ii) dividing the measurement spectral data set D'^jiio)) amplitude at a particular 
frequency in the measurement spectral data set Di(a)) , 

iii) by the amplitude of the reference measurement [(o) at the same frequency, 

b) to create a normalized data wavefield XjJ (o?) . 

12. The method of claim 4 wherein the error bound of the minimizing step is calculated 
by first assigning an initial error E^, then dividing each subsequent iterated error 

bound by the initial error . 

13. An iterative method of using normalized waveform inversion to obtain a model 
parameter describing one or more physical properties of a medium, the method 
comprising the steps of: 

a) inputting one or more each of NS source and NG measurement spatial locations; 
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} 

b) measuring time domain data at each of the NG measurement locations resulting 
from an input waveform at one of the NS source locations propagating through a 
medium, 

i) for each of the NS source locations, 

ii) thereby forming a time domain measurement data set 

Z>J (0, (/ = 1 NG; i = NT) , where d denotes a data tensor; 

c) Fourier transforming the time domain measurement data set />^ (^) to create a 
measurement spectral data set D'^jiiw) having frequency and amplitude 
information for each of the NG measurement locations; 

d) normalizing the measurement spectral data set Di{co) to create a normalized data 

wavefield Tj^io)); 

e) modeling the medium using an iterated model parameter m describing one or 
more physical properties of the medium, 

the A^5 source and the NG measurement spatial locations used as respective 
model input and model response spatial locations contained within the model 
of the medium, 

and initializing the iterated model parameter m with corresponding one or 
more known bulk properties of the medium being modeled, 
said modeling step comprising: 

i) creating a measurement model by: 

(1) applying a delta function source collocated with the NS source, 

(2) modeling the response at the NG measurement locations, using a time 
domain modeling method, to create a synthetic medium response at the 
receiver due to the source, P^7 (0 ' denoting the model response, 

(3) repeating the applying and modeling steps at each of the NS source 
locations and NG measurement locations until the measurement model is 
full, and 

ii) Fourier transforming the model response (^) , to obtain a frequency 
domain synthetic response P^7(^) ^ 
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iii) forming a normalized modeled wavefield using the frequency domain 
synthetic response T^lio)) = P;(ty)[Pj7(iy)]"' ; and 

f) minimizing a weighted error, between the normalized data wavefield T^^(^) and 
the normalized modeled wavefield T]l(co) of the response of the medium, to a 

level below an error bound, 

g) said weighted error met by using the iterated model parameter m , known as the a 
minimized model parameter m . 

14. The method of claim 13 further comprising the step of; 

a) storing as a computer program in at least one computer-readable medium, 
i) said steps of claim 13. 

15. An article of manufacture for normalized waveform inversion using a computer, said 
article comprising: 

a computer readable medium comprising instructions for a computer to execute, said 
execution comprising the steps of: 

a) inputting one or more each of NS source and NG measurement spatial locations; 

b) measuring time domain data at each of the NG measurement locations resulting 
from an input waveform at one of the NS source locations propagating through a 
medium, 

i) for each of the NS source locations, 

ii) thereby forming a time domain measurement data set 

D^ji (0» (/ = 1 ^G; / = 1 - , where d denotes a data tensor; 

c) Fourier transforming the time domain measurement data set D^^(t) to create a 
measurement spectral data set D^,.(ft;) having frequency and amplitude 
information for each of the NG measurement locations; 

d) normalizing the measurement spectral data set D^^.(£2;) to create a normalized data 

wavefield T^((y); 
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e) 



f) 



modeling the medium using an iterated model parameter m describing one or 

more physical properties of the medium, 

the NS source and the NG measurement spatial locations used as respective 
model input and model response spatial locations contained within the model 
of the medium, 

and initializing the iterated model parameter m with corresponding one or 
more known bulk properties of the medium being modeled, 
said modeling step comprising: 

i) creating a measurement model by: 

(1) applying a delta function source collocated with the NS source, 

(2) modeling the response at the NG measurement locations, using a time 
domain modeling method, to create a synthetic medium response at the 
receiver due to the source, PJ {t) , m denoting the model response, 

(3) repeating the applying and modeling steps at each of the NS source 
locations and NG measurement locations until the measurement model is 
full, and 

ii) Fourier transforming the model response (r) , to obtain a frequency 
domain synthetic response V]l(a)) ; 

iii) forming a normalized modeled wavefield using the frequency domain 



minimizing a weighted error, between the normalized data wavefield Tj-ico) and 
the normalized modeled wavefield T^iio)) of the response of the medium, to a 
level below an error bound, 

i) said weighted error met by using the iterated model parameter m , known as 
the a minimized model parameter m . 

iterative method of using normalized waveform inversion to obtain a model 
ameter describing one or more physical properties of a medium, the method 
iprising the steps of: 
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a) inputting one or more each of NS source and NG measurement spatial locations; 

b) measuring time domain data at each of the NG measurement locations resulting 
from an input waveform at one of the NS source locations propagating through a 
medium, 

i) for each of the NS source locations, 

ii) thereby forming a time domain measurement data set 

D^ji U)* 0 = 1'" / = 1 - iV5) , where d denotes a data tensor; 

c) Fourier transforming the time domain measurement data set D^ (0 to create a 
measurement spectral data set D^,.(fi£;) having frequency and amplitude 
information for each of the NG measurement locations; 

d) normalizing the measurement spectral data set D^io)) to create a normalized data 

wavefield Tji(a)); 

e) modeling the medium using an iterated model parameter m describing one or 
more physical properties of the medium, 

the NS source and the NG measurement spatial locations used as respective 
model input and model response spatial locations contained within the model 
of the medium, 

and initializing the iterated model parameter m with corresponding one or 
more known bulk properties of the medium being modeled, 

said modeling step comprising: 

i) creating a measurement model by: 

(1) applying a delta function source collocated with the source, 

(2) modeling the response at the NG measurement locations, using a 
frequency domain modeling method, to create a frequency domain 
synthetic response at the receiver due to the source, P^^ (^) » ^ 

denoting the model response, 

(3) repeating the applying and modeling steps at each of the NS source 
locations and NG measurement locations until the frequency domain 
synthetic response is full, and 
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ii) forming a normalized modeled wavefield using the frequency domain 
synthetic response T;(iy) = P;((y)[P,7((y)]"* ; and 

f) minimizing a weighted error, between the normalized data wavefield TjJ (cd) and 
the normalized modeled wavefield TJI(a)) of the response of the medium, to a 
level below an error bound, 

i) said weighted error met by using the iterated model parameter m , known as 
the a minimized model parameter m . 

17. The method of claim 16 further comprising the step of; 

a) storing as a computer program in at least one computer-readable medium, said 
steps of claim 16. 

18. An article of manufacture for normalized waveform inversion using a computer, said 
article comprising: a computer readable medium comprising instructions for a 
computer to execute, said execution comprising the steps of claim 16. 
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